Mon. Not. R. Astron. Soc. 000, [THIS] (2008) Printed 29 January 2010 (MN WI^ style file v2.2) 



A near- infrared study of AGB and red giant stars in the 
Leo I dSph galaxy 

E. V. Held^*, M. Gullieuszik\ L. Rizzi^, L. Girardi\ P. Marigo^ and I. Saviane^ 

^ INAF/Osservatorio Astronomico di Padova, vicolo dell'Osservatorio 5, 1-35122 Padova, Italy 
^ Joint Astronomy Centre, 660 N. A'ohoku Place, University Park, Hilo, HI 96720, USA 
^ Dipartimento di Astronomia, Universita di Padova, vicolo dell'Osservatorio 2, 1-35122 Padova, Italy 
* European Southern Observatory, Casilla 19001, Santiago 19, Chile 

^ : 

o\ '. ^ 

• Accepted . . . Received . . . ; in original form . . . 

ABSTRACT 

A near-infrared imaging study of the evolved stellar populations in the dwarf 
spheroidal galaxy Leo I is presented. Based on JHK observations obtained with the 
WFCAM wide-field array at the UKIRT telescope, we build a near-infrared photo- 
metric catalogue of red giant branch (RGB) and asymptotic giant branch (AGB) 
stars in Leol over a 13.'5 square area. The V — K colours of RGB stars, obtained by 
combining the new data with existing optical observations, allow us to derive a dis- 
tribution of global metallicity [M/H] with average [M/H] = —1.51 (uncorrected) or 
[M/H] = — 1.24±0.05 (int)±0.15 (syst) after correction for the mean age of Leol stars. 
This is consistent with the results from spectroscopy once stellar ages are taken into 
account. Using a near- infrared two-colour diagram, we discriminate between carbon- 
and oxygen-rich AGB stars and obtain a clean separation from Milky Way foreground 
stars. We reveal a concentration of C-type AGB stars relative to the red giant stars in 
the inner region of the galaxy, which implies a radial gradient in the intermediate-age 
(1-3 Gyr) stellar populations. The numbers and luminosities of the observed carbon- 
and oxygen-rich AGB stars are compared with those predicted by evolutionary models 
including the thermally-pulsing AGB phase, to provide new constraints to the models 
for low-metallicity stars. We find an excess in the predicted number of C stars fainter 
than the RGB tip, associated to a paucity of brighter ones. The number of 0-rich 
AGB stars is roughly consistent with the models, yet their predicted luminosity func- 
tion is extended to brighter luminosity. Although these discrepancies can be partly 
ascribed to significant uncertainties in the Leol star- formation history and incom- 
pleteness of the spectroscopic samples of C stars fainter than the RGB tip, it appears 
more likely that the adopted evolutionary models overestimate the C star lifetime and 
underestimate their /•sT-band luminosity. 

Key words: Galaxies: individual: Leol - Galaxies: stellar content - stars: AGB and 
post-AGB - stars: carbon - Local Group 



1 INTRODUCTION 

The dwarf spheroidal galaxy Leol was the first shown 
to be dominate d by a relatively young stellar population 
jLee et al.ll 19931 '). HST deep photometry confirmed the pres- 
ence of a pr ominent intermediate-age stellar popula tion and 
young stars (|Gallart et al. l ll999l : ICaputo et al.|[l999l ). Quan- 
titative analysis of the HST data established that most star 
formation activity (70%-80%) occurred between 7 and 1 
Gyr ago, with some residual star formation until at least 
~ 300 Myr ago. The detection of a blue horizontal branch 
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jHeld et al.lbOOOl ') provided the first evidence for an old stel- 
lar population, later confirmed by the discovery of RR Lyrae 
variable stars (Held et al. 200^). 

The presence of an intermediate age population in Leo I 
was also revealed by an extended asymptotic giant branch 
(AGB). After the spectro s copic discovery of a carbon star in 
Leol by Aaronson et all (Il983|). JHK follow-up photome- 
try by lAaronson fc Mouldl ( 1985h showed that the carbon 



stars (4 newly identified on the basis of their red near- 
infrared colours) belong to an extended AGB typical of an 
intermediate-age stellar population. The presence of a sig- 
nificant population of C stars (19 stars) in Leol was con- 
firmed by objective grating spectroscopy in the blue-green 
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spect ral region (Swan C2 bands) by I Azzopardi et alj (| 19851 . 
1 19861 ). Using a 4- filter technique based on the RI filters and 
CN an d TiO intermediate-band filters. [Pemers fc BattineilH 
further surveyed Leo I for C stars, i dentifying 13 C- 
star candidates, 6 of which in common with I Azzopardi et al.l 
(| 19861 ). iMenzies et al] (|2002l ) obtained near-infrared (NIR) 
photometry to Ka = 16 in a 7'.2 square field allowing a study 
of the AGB content of Leo L They showed that most lumi- 
nous AGB stars are carbon stars, several of which are vari- 
able, and detected the presence of 5 very red stars near the 
AGB tip, probably young (~ 2 Gyr old) dust-enshrouded 
stars, also candidates for Mira variability. 

Several spectroscopic studies recently a ddressed the 
metal licity and chemical evolution of LeoL iBosler et al.l 
used the Can triplet method to analyse Keck-LRIS 
spectra of 102 red giant branch (RGB) stars, and found 
[Fe/H] = -1 34. A mean meta llicity [Fe/H] = -1.31 was 
measured by iKoch et al.l (|2007l ) for 58 red giants. In both 
studies, the metallicity distribution is well described by a 
Gaussian function with a la width of 0. 25 dex, and a full 
range in [Fe/H] of approximately 1 dex. In lGuUieuszik et al.l 
(|2009D we obtained metallicities for 54 stars in Leo I from 
measurements of the Ca triplet line strengths, and found a 
mean metallicity [Fe/H] = -1.4 and [M/H] = -1.2 on a 
scale tied to the [M/H] ranking of Galactic globular clus- 
ters (GGCs). A narrow and symmetric metallicity distri- 
bution, with a very low instrinsic metallicity dispersion, 
c"[M/H] = 0.08 dex, appears consistent with a prompt initial 
enrichment and galactic winds expelling the metals products 
of stellar evolution. We also suggested the presence of a ra- 
dial age gradient, with RGB stars in the inner part of Leo I 
being on average younger than those in the outer regions. 

In this study, we further investigate the star-formation 
history of Leo I and the presence of stellar population gra- 
dients using deep NIR wide-field observations and analysing 
the AGB stars as tracers of intermediate-age stellar popu- 
lations. Our NIR observations provide new constraints on 
the star-formation history of the galaxy, which are set by 
the number, type (carbon vs. ox ygen-ri ch), and luminos i- 
ties of AGB stars (see, e.g.. lCioni 2007: Groenewegenll2007^ . 
Near-infrared photometry is ideal to this purpose because 
the cool AGB stars emit most of their flux at near-infrared 
wavelengths. The most evolved, dust-enshrouded AGB stars 
can be detected only at infrared wavelengths. In addition, 
the NIR photometry of red giant stars can be used to provide 
independent information on basic stellar properties. 

Together with a companion paper on Leo II 
iGuUieuszik et al] I2OO8I ). this study is part of a survey 
of evolved stellar populations in dwarf galaxies in the Local 
Group and its vicinity, aimed at providing an observational 
database for the calibration of AGB stars models in a wide 
range of metallicities as a tool for a reliable interpretation 
of observations. For distant resolved systems (e.g., Virgo 
cluster galaxies observed with JWST or ELT) this will 
represent the main observational fact for modelling the 
galaxy's star formation rate in the last few Gyr. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of Leo I were carried out on A pril 19-20, 2005 
using the NIR wide field camera WFCAM (|Casah fc et all 



Table 1. Observing log. 



Filter 


N- 

' ima 


DIT(s) 


N 

^ ' exp 




Microsteps 


J 


6 


5.0 


2 


9 


2x2 


H 


6 


5.0 


2 


9 


2x2 


K 


10 


5.0 


2 


9 


2x2 



Table 2. The NIR catalogue of Leo I stars over WFCAM ar- 
ray No. 2. A few lines are shown here for guidance regarding its 
form and content, while the full catalogue is available from the 
electronic edition of the journal. 



ID 


a 


(J2000) 


<5 (J2000) 


J 


H 


K 


14476 


10 


08 


29.36 


-(-12:12;05.6 


19.48 


18.60 


18.43 


14447 


10 


08 


03.90 


-1-12:12:06.8 


19.64 


18.93 


18.93 


200324 


10 


08 


54.77 


+12:12:06.9 


18.78 


18.26 


18.32 


14467 


10 


08 


41.75 


-1-12:12:07.3 


20.78 


19.95 


20.05 


22945 


10 


08 


38.14 


+12:12:07.7 


20.64 


19.85 


20.73 



l2007h at the UKIRT at Mauna Kea, Hawai'i. WFCAM is 
equipped with 4 Rockwell Hawaii-II HgCdTe detectors with 
a scale of 0'.'4 pixel" \ The size of each detector is 2048 x 2048 
pixels, corresponding to 13'6 x 13'6, with a space between 
the arrays of 94% of the detector size. Sirice the tidal radius 
of Leo I is 12!6 (|lrwin fc Hatzidimitrioul Il995h . one array 
(No. 2) was sufficient to study Leo I. The remaining detec- 
tors provide a useful estimate of the background/foreground 
field contamination. 

The observing strategy and ins trumental setup are th e 
same as used for our study of Leo II (|Gullieuszik et al]|2008l ) . 
The observations, summarised in Table [TJ consisted of 4 
micro-stepped 10 s images (two 5 s coadds) on a 9-points 
jitter pattern, yielding a 360 s exposure time. In total, 6 
observations were obtained in J and H and 10 in the K 
band, yielding on-target integration times of 36min in J 
and H and 60 min in K. The 2x2 "small" microstepping 
was used to obtain a better spatial sampling. The raw data 
were calibrated using the WFCAM p ipeline provided by the 
VISTA Data Flow System Project (|Dve et al.ll2006h . The 
product of the pipeline are 4kx4k oversampled combined 
images ( "leavstacks" ) with a spatial resolution 0'.'2 pixel" ^. 

Our analysis was limited to the array centred on Leo I 
(No. 2) and a second array (No. 3) used as a control field 
to estimate the contamination by foreground Galactic stars 
and background galaxies. Object detection and point-spread 
function (PSF) photometry was performed o n each "leavs- 
tack" image using daophot and allframe (|Stetson|[l987l . 
Il994h . The PSF was generated with a "Penny" function with 
quadratic dependence on the position on the frame. The as- 
trometric calibration was performed by the pipeline using 
the 2MASS Point Source Catalogue (|Skrutskie et al.ll2006l . 
PSC) as a reference, with a final absolute systematic accu- 
racy of the order O'.'l. The positions of the sources in the 
final JHK catalogue were converted from pixel units to the 
J2000 equatorial system using this astrometric calibration 
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Fig ure 1. NIR and opti cal-NIR colour-magnitude diagrams of Led. Superimposed on the {J — K), K CMD are theoretical isochrones 
fro mlGirardi et al.l 1 I2OO2I) . The age is 6 Gyr and the metall icities bracket that of Leo I: Z = 0.002 and Z = 0.0004. Carbon stars identified 
bv lAzzopardi et aLlTTosg) andlDemers fc Battinelll 1 I2OO2I) are shown as diamonds and squares, respectively. The starred symbols arc the 
very red stars identified bv lMenzies et ah I l l2002t) . 



and an updated version of the IRAF0 package that included 
support for the ZPN projection adopted by the pipehne. 

Our raw photometric catalogue was transformed to the 
2MASS photometric system by applying the colour terms 
between WFCA M and 2MASS systems derived for the 
UKIDSS surveys (|Dve et al.11200^ . and zero points directly 
computed from the median of the JHK magnitude differ- 
ences between 2MASS and our WFCAM photometry in the 
colour range 0.0 <J—K< 1.0, for stars in common with the 
2MASS/PSC. Only stars having 2MASS S/N ratio > 10 and 
photometric errors < 0.1 mag were used, which essentially 
limits the comparison to stars brighter than the RGB tip. 
The r.m.s. of the residuals is 0.04 mag in the J band and 0.05 
mag in the H and K bands. We verified that the zero points 
would not c hange significan t ly (> 0.01 mag) using the colour 
terms from iHodgkin et al.l (|2009l ). We note that the colour 
terms were defined in the colour range 0.0 <J—K< 1.0. 
Outside this interval (e.g., for very red stars), our measure- 
ments should be considered on a (linearly transformed) WF- 
CAM system rather than on the 2MASS system. The final, 
calibrated NIR photometric catalogue, which includes star- 
like objects detected in at least 2 images in 2 bands, is pub- 
lished in the electronic issue of the journal. A few lines are 
presented in Table [2] to illustrate its content. 

For the completeness and photometric errors we as- 
sumed th«_j;esults_o£Mtifi£i^ stars experiments used for 
Leo II jGullieuszik et al]|2008l ) , since the observing strategy, 

^ The Image Reduction and Analysis Facility (IRAF) software 
is provided by the National Optical Astronomy Observatories 
(NOAO), which is operated by the Association of Universities 
for Research in Astronomy (AURA), Inc., under contract to the 
National Science Foundation. 



exposure times, and sky conditions were identical for the two 
galaxies, and the crowding is very similar and low in the 
NIR. The completeness factor is > 50% for stars brighter 
than if ~ 20. All results of our analysis are based on pho- 
tometry of stars brighter than K — 18, for which we have 
completeness factor ~ 100% and internal photometric errors 
< 0.02 mag. 



3 COLOUR-MAGNITUDE DIAGRAMS 

Figure [1] presents our (J — K), K colour-magnitude dia- 
gram (CMD) of the Leo I field, along with a diagram com- 
bining NIR data and optical photometry obtained with 
the EMMI came ra at the ESO/La Silla NTT telescope 
iHeld et al.ll2000l ). Our photometric catalogue was cleaned 
using the sharp parameter calculated by ALLFRAME to re- 
move noise peaks, diffuse objects, and other spurious detec- 
ti ons. A magnitude - depend ent sharp cutoff was adopted, as 
in lCullieuszik et all (j2007bl . their figure 2). 

Our dataset provides a nearly complete sample of 
upper-AGB stars belonging to the intermediate-age stellar 
populations of Leo I, with a red tail of C stars which e xtend s 
to (J— A') ~ 3.5, in agreement with iMenzies et al.l (|2002h . 
In addition, the NIR CMD of Leo I samples the RGB from 
the tip (clearly visible at K ^ 16) down to K ^ 2 0.5. The 
RGB is comprised between the lGirardi et all (|2002l ) theoret- 
ical isochrones for an age 6 Gyr and metallicities Z — 0.002 
and Z = 0.0004. 

The cool AGB C-stars (mostly variable) are the bright- 
est objects in Leo I in the K band. At optical wavelengths, 
their luminosity dramatically decreases and the reddest 
among them are generally undetected in the B band. The 
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Figure 2. A 1' X 1' region in Leo I as seen in the V band by (NTT/EMMI data, left panel) and in the K band (WFCAM data, right 
panel). Open circles mark stars brighter t han K = 18, while o pen squares show the position of the stars "A" (on the left side of the 
image) and "B" (on the right) identified bv lMenzies et ah These stars are very bright in the near infrared and yet undetected at 

optical wavelengths. 



completeness of the optical catalogues is also complicated by 
the large amplitude of luminosity variations in the optical 
domain. Still, many upper- AGB stars are present in the op- 
tical catalogues, although some are very faint (see Fig.[T]). As 
a result, archival optical catalogues may be used to retrieve 
AGB stars detected in the NIR (or mid-infrared) and char- 
acterise their SEDs. We cross-identified the stars in our NIR 
catalo gue with the carbon star s urvey s of lAzzopardi et al] 
l|l98d) and lDemers fc BattineUH (120021) a nd with the list of 
very red objects in iMenzies et al! ( 20021 ) . using positional 
coincidence within a 2'.'5 radius. Star C IO was incorrectly 
ide ntified by Demers fc BattineliU (|2002l ) with star no. 7 in 
the lAzzopardi et al. ( 1986h catalogue. Since the distance be- 
tween D02-C10 and ALW-7 is 32" , and both stars are located 
within 2'.' 5 from two distinct sources in our catalogue (with 
colours and magnitudes typical of C stars), we conclude that 
they are distinct sta rs. All of the very red stars detected by 
iMenzies et al] (|2002l ) within our field-of-view (namely, stars 
A, B, C, and D) are also present in our NIR CMD (Fig.Q. 
Among them, only star D is visible in the optical bands. 
As an example. Fig. [2] shows two of the reddest stars. They 
are very bright in the K band, yet fall below the detection 
threshold in the V image. 

Different techniques adopted for s earching C stars lead 
to sa mples with different biases (see iBattinelli fc DemersI 
I2OO9D . The reddest stars are too faint in the optical spec- 
tral range to be detected with objective grating tech- 
niques, while the 4-filter t ec hnique misses th e blue st stars 
l|Demers fc Battinellil [2OO2I ). lAzzopardi et al.l (|l986l ). using 
an objective grating operating in the 5000 — 7000 A range, 
did not detect the brig htest, reddest C stars (as noticed 
by IMenzies et al] 12003 ) . while iDemers fc Battineilil (|2002l ) 
could not retrieve the bluest, i.e. warmest, spectroscopically 
identified C stars. 



In Sect. |6] we will present a method for selecting C- and 
O-rich AGB stars based on two-colour NIR diagrams. Our 
selection will include all previously known C-stars and, in 
addition, will allows us to separate the foreground Milky 
Way stellar population. 



4 DISTANCE 

The luminosity of the RGB tip (TRGB) in the / band 
has long been used as a valuable distan ce indicator 
(|Da Costa fc Armandrogil990l : lLe"e et al.lll993l ). The depen- 
dence of the TRGB luminos ity on age and metallicit y is min- 
imised in the 7-band (e.g., ISalaris fc Girardill2005h . In the 
NIR, the luminosity of the tip depends on age and metallic- 
ity in a more complex way. For instance, the if-band magni- 
tude of intermediate-age stars at the TRGB is fainter than 
that of old stars; while the TRGB K luminosity ris es with 
increasing metallicity (e.g., ISalaris fc GirardU |2005^ . For a 
population that becomes more metal-rich with time as a re- 
sult of galaxy chemical evolution, the two effects may p artly 
balance. As discussed in lGuUieuszik et al.l (|2007al . [2008h . dis- 
tance estimates consistent with those obtained from opti- 
cal methods can be derived from NIR photometry, provided 
that the star-formation history (SFH) of the galaxy is at 
least approximately known. Testing our ability to measure 
distances from NIR data is important in view of the future 
generation of ground based and space telescopes that will op- 
erate mainly or exclusively at infrared wavelengths (JWST, 
adaptive optics at the ELT). 

In a complex stellar population, the RGB is generally 
wide and the cutoff does not have a con stant luminosity 
(see, e.g., the discussion of Fornax dSph bv lWhitelock et"al] 
[iooi). However, the TRGB can still be operationally de- 
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Table 3. Observed magnitude of the TRGB and corrected dis- 
tance modulus fiQ derived for Leo I from JHK photometry. The 
distance was population-corrected by comparing the absolute 
magnitudes of the TRGB calculated from synthetic CMDs for 
the population mix of Leo I (Af^") and a pure old population 
(M°'<^) (see text for details). 



band 


^TRGB 






MO 




J 


17.00 


-5.310 


-5.326 


22.09 


0.18 


H 


16.30 


-6.179 


-6.329 


22.01 


0.19 


K 


16.14 


-6.293 


-6.340 


22.01 


0.20 



fined to represent the mix population, and compared with 
synthetic CMDs that refiect the SFH of the system. The 
method is employed here to measure the distance to Leo I. 
An objective estimate of the magnitude of the TRGB was 
obtained by fitting the A'-band luminosity function with a 
step function convolved with a Gaussian kernel represen- 
tative of the photometric errors. This method, extensively 
applied by our group (e.g.. iMomanv et al.ll2002 ). has given 
in the I band results consisten t within la with the Max- 
imum Likelihood Algorithm of iMakarov et al.l (|2006l ) (see 
iRizzi et al.ll2007l '). 

The magnitudes of the TRGB measured in the JHK 
bands are listed in Table O The errors associated to these 
magnitudes are essentially the uncertainty on the absolute 
photometric calibration, since the uncertainties due to the 
fitting algorithm are less than 0.01 mag and the internal 
photometric errors at the level of the TRGB are negli- 
gible (see Sect. [2]). To derive the distance to Leo I, the 
JHK TRGB magnitudes were compared with the empir- 
ical calibrations of M™'^^ as a function of [M/H] from 
IValenti et all (|2004bh . This calibration, whose intrinsic sys- 
tematic error is (t = 0.16 mag, is based on Galactic globu- 
lar clusters, i.e. old stellar populations. Applying it to the 
complex stellar populations of Leo I requires a "population 
correction" based on synthetic CMDs to take into account 
its star formation history. The distance to Leo I was then 
computed as 



(m - M)o ^mx-Ax- (Ma + AA/^) 



(1) 



where mx is the observed TRGB magnitude. Ax is the ex- 
tinct i£n^_JVfA_Js_Uie absolute TRGB magnitude as derived 
from lValenti et al] (|2004bl ) calibration, and AMx is the pop- 
ulation correction needed to account for the effects of the 
different SFH of Leo I and the GGCs. 

The population correction AMa was calculated using a 
simulated CMD of Leo I with constant metallicity to build 
the JHK luminosity function (LF) and measure the theo- 
retical TRGB magnitude, M^" (see lCuUieuszik et al.|[2007al 
for details). Our measurements were then repeated by se- 
lecting only stars older than 10 Gyr (M^''*). The population 
correction was then AMa = M^" - Mf"^. The values are 
listed in Table [3] The population effect is positive for all 
bands, i.e. the RGB tip for the overall stellar population is 
slightly fainter than that of an old population. 

Our distance estimates for Leo I, computed from Eq. ([l]) 
with a reddening Eb-v ~ 0.03, are given in Table [31 The 
weighted mean of the distances derived from the JHK 
bands is 22.04 mag (corresponding to 256 kpc) with a stan- 




(V-K)o 

Figure 3. Derivation of a photometric Metallicity Distribution 
Function of RGB stars in Leo I using V — K colours. Upper panel: 
the metallicity distribution of red giants in Leo I. The solid line 
is a fit of a Gaussian function to the data with [M/H] > —1.8. 
The method is illustrated in the lower panel: the stars on the 
upper RGB of Fornax are interpolated across analytical fits ( dot- 
ted lines) to the RGB fiducial lines of template Galactic globular 
clusters. The thick line marks the best-fit mean metallicity. 



dard deviation of 0.11 mag. The errors were calculated by 
error propagation from the photometric calibration uncer- 
tai nty (see Sect. and the intrinsic systematic error on 
the lValenti et all ([2004bl 'l calibration (0.16 mag) This esti- 
mate agrees with the results of optical studies. iHeld et"al] 
(j200lh found (m - M)o = 22. 04 ± 0.14 from the V m agni- 
tude of Leo I RR-Lyrae stars; iBellazzini et al.l (|2004h mea- 
sured (m — M)o = 22.02 ±0.13 f rom the I magnitude of the 
TRGB. Using the same method. [Mendez et alj|j2002h found 
(m-M)o = 22.05 ±0.10 (internal) ± 0.18 (systematic). Our 
result therefore confirms the possibility of reliably estimat- 
ing the distance to resolved stellar systems with complex 
stellar populations using NIR photometry alone. 



5 METALLICITY 

5.1 Metallicity distribution 

The metallicity distribution function (MDF) of RGB stars 
is an important constraint for models of chemical evolution 
of dwarf galaxies. We derived a photometric MDF for RGB 
stars in Leo I by comparing their V — K with the colour of 
GGCs of known met allicity, following a te chnique developed 
for optical colours bv lSaviane et al.l (|2000l ). Stars on the up- 
per RGB were interpolated in the (V — K), Mk diagram 
across analytical fits to the RGB fiducial lines of template 
Galactic globular clusters. We assumed for Leo I a distance 
of (m - M)o = 22.04 and a reddening Eb-v = 0.03. In 
our implementation of the method, we employed [M/H] val- 
ues for a set of GGCs spanning a wide range in metallic- 
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Figure 4. Comparison of the MDFs of Leo I RGB stars 
derived from NIR pho tometry (this work) and spectroscopy 
llGuUieuszik et al 1 l200g|) . The photometric MDF was corrected 
for age effects assuming a mean age of 5 Gyr. 



ity, from IValenti et al] (l2004al'l. The complete list of GGCs 
with metallicities is given by Gullieuszik et all (|2007ah . The 
global metallicity [M/H], which measures the abundance of 
all heavy elements, is the most appropriate to estimate the 
metallicities of dwarf spheroidal galaxies (having [a/Fe] ra- 
tios close to solar) by comparison with the photometric prop- 
erties of Milky Way globular clusters, which generally show 
an overabundance of a-elements relati ve to iron that is a 
function of the cluster metallicity (see ICeisler et alj |2007| . 
and refs. therein). 

The resulting photometric MDF of Leo I, obtained 
for red giant stars down to ~ 1 mag below the TRGB 
(—6.0 < Mk < —5.0), is presented in Fig. [3] (upper panel). 
With our choice, the MDF shows negligible contamination 
by stars on the blue side of the RGB, which probably in- 
clude both younger stars in Leo I and Galactic foreground. 
The distribution has an average [M/H] = —1.51 and is quite 
well described by a Gaussian centred at [M/H] = —1.45 with 
a measured dispersion 0.19 dex. The standard deviation of 
the mean, 0.05 dex, is adopted as internal error of the mean 
metallicity. 

Since the stellar populations in Leo I are on average 
younger that the Galactic globular clusters, we computed a 
correction to the mean metallicity of Leo I by adopting the 
simpli fied approach that a ll RGB stars have an age of 5 Gyr 
(e.g., iGallart et al]|l999l '). Given two RGB stars with the 
same metallicit y, the younger has b luer colours. By interpo- 
lating over the lGirardi et al.l (|2002l ) isochrone set, we found 
that at Mk ~ —5.0 a star of 5 Gyr is more metal rich by 
0.27 dex than a 12.5 Gyr old star having [M/H] = -1.51 
and the same V — K colour. By applying this correction, 
the mean metallicity of Leo I from NIR photometry is then 
[M/H] = -1.24 ± 0.05 (random) ±0.15 (systematic). The 
systematic error was estimated from the uncertainty on the 



Figure 5. Individual metallicity measures as derived from spec- 
troscopy and from photometry, without any age correction. Dif- 
fer ent symbols are used fo r stars with different ages, as derived 
by iGuUieuszik et aL The solid line is the bisector, while 

the dashed line corresponds to the age correction for 5 Gyr old 
stars. The error bars represent the mean error of the spectro- 
scopic metallicities and the r.m.s. (internal) error of [M/H] from 
optical-NIR colours. 



photometric zero point of our calibration, which is 0.07 mag 
on the V — K colour (the quadratic sum of the V and K 
zero-point uncertainties). Shifting the Leo I RGB by ±0.07 
mag in colour results in a ±0.15 dex variation in metallicity. 

The mean metallicity of red giant stars in Leo I 
in our spectroscopic an alysis is [M/H] = —1.26 ± 0.16 



Gullieuszik et all I200S ) . Prev ious work by iBosler et al.l 
20071 ) and iKoch et all (|2007l ) ed to mean metallicities 
[Fe/H] = -1.34 and [Fe/H] = -1.31, respectively. Our pho- 
tometric age-corrected metallicity therefore agrees with all 
recent spectroscopic measures. 



5.2 Comparison with spectroscopy 

The metallicities of individual RGB stars in Leo I, derived 
from NIR photometry, can be directly compared with the 
results of spectroscopic measurements. In Fig. [4] the metal- 
licity distribution of Leo I RGB stars inferred from V — K 
colours is com pared with the spectros copic metallicity dis- 
tribution from I Gullieuszik et al.l (|2009l ) (we refer the reader 
to that paper for a full account of previous spectroscopic 
work). The photometric MDF was approximately corrected 
for age population effects, assuming a common age of 5 Gyr 
for Leo I stars. The [M/H] distributions are similar, both in 
the mean values and the range. The primary difference is in 
the low-metallicity tail of the photometric MDF. This can 
be partially explained by a fraction of contaminating AGB 
stars, although the main reason is probably the presence 
of young stars with colour bluer than the RGB, which are 
incorrectly treated as metal-poor stars. 
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Figure 6. Near-IR two-colour diagrams of stars brighter than the RGB tip in the Leo I field {left panel) and tiic control field (right 
panel), with superimposed the regions used to distinguish stars in Leo I from foreground stars. Region (1) delimits probable O-rich AGB 
stars in Leo I, regions (2) and (3) contain Galactic dwarf stars, region (4) is populated by C stars in Leo I. Different symbols and colours 
are used for stars in different regions. Also shown are the loci of giant stars [dashed line) and main-sequence stars (solid line) from 
iBessell fc Brett! lll988l) . Note the absence of stars belonging to Leof in the control field, which is located beyond the "tidal radius" of 
Leo I. 



The significant over lap between our sampl e of RGB 
stars and the dataset of iGullieuszik et alj (|200y) (39 stars 
in common) also allows a direct comparison of photometric 
and spectroscopic metallicity estimates on a star-by-star ba- 
sis. Individual photometric and spectroscopic measures are 
plotted in Fi g. [5l where stars in diff erent age intervals (as 
estimated by IGullieuszik et aill2009l) . are marked with dif- 
ferent symbols. In this plot, the metallicities derived from 
photometry have not been corrected for the age difference. 
Were all Leo I stars as old as the GGCs, we would expect 
the spectroscopic and photometric metallicities to be equal 
within the errors, and cluster about the bisector in Fig. (5) 
This is clearly not the case, as most stars are on the left of 
the unity relation, i.e. their metallicities are underestimated 
by the V—K colours. This is consistent with the mean age of 
the bulk of the Leo I stellar populations. We have plotted an 
age-corrected relation (dashed line in Fig. [5} , by calculating 
the metallicity correction to be applied to our photomet- 
ric measures for young stellar populations (A[M/H] — 0.27 
dex, as above). Indeed, many stars are close to this line, 
with the youngest stars located even more to the left, as 
expected. The observed shift inV — Koi young stars is con- 
sistent with our age determination based on optical colours. 
The metallicities derived using the V — K colours agree with 
the spectroscopic estimates, once stars' ages are accounted 
for. Alternatively, the V — K colours could be used to es- 
timate the ages of stars once their metallicities are known 
from spectroscopy. 



6 AGB STARS 

We now turn our attention to the AGB star content of Leo I. 
Our WFCAM JHK catalogue provides the database for 
a study of the intermediate-age AGB stars in Leo I with 
unprecedented photometric accuracy and spatial coverage. 
From this database, we now select new samples of C-rich 
and O-rich AGB stars and discuss their properties. 



6.1 Two-colour diagrams 

The two-colour d iagram is an efficient tool to separate dwarf 
and giant stars (|Aaronson fc Mouldl [l985l : iBessell fc Bret^ 
Il988h . As such, it is particularly useful when studying ex- 
tragalactic systems for separating the stellar populations of 
the system from the Galactic foreground. 

Figure [S] shows the NIR two-colour diagrams of stars 
brighter than K = 16.0 in the field of Leo I and in a compar- 
ison field beyond the "tidal radius" of the galaxy. The mag- 
nitude cutoff was chosen to include only stars brighter than 
the RGB tip, leaving out first-ascent red giant stars. Fol- 
lowin g the approach of our Leo 11 study jGullieuszik et ahl 
|2008| ) , we defined 4 regions in the two-colour diagram aimed 
at discriminating C stars from O-rich AGB stars in Leo I, 
both against the Galactic foreground. The star counts in 
the different regions are given in Fig. |S] 

Regions 2 and 3 contain the loci of gian t and main- 
sequence Galactic stars from iBesseU fc Bret^ (|l988l ). The 
star counts are equal in Leo I and the control field within 
the statistical y/N fluctuations, indicating a population of 
foreground Milky Way stars. Regions 1 and 4 are popu- 
lated only in the case of the Leo I field, which proves that 
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Figure 7. The position of AGB stars and foreground stars, clas- 
sified according to tfie NIR two-colour diagram of Fig. [6] in the 
CMD of the Leo I field. O-rich AGB stars (region 1, red squares) 
and C stars (region 4, yellow filled circles) in Leo I are effectively 
discerned from foreground Milky Way stars (region 2, green open 
squares and region 3, cyan open circles). The solid curve is a 
mean col our-magnitude rela tion for C star in Local Group dwarf 
galaxies llTotten et al.ll2000l) . scaled to the distance of Loo I. 



stars in those regions belong to the dSph galaxy. The stars 
in region 1, which are close to the iBessell fc BrettI (|l988l ') 
line for giant stars, are identified with O-rich AGB stars in 
Leo I. The photometry of candidate O-rich AGB stars in re- 
gion 1 is given in Table 3] Region 4 contains stars which 
become increasingly reddened in both J — H and H — K 
colours. In dSph g alaxies this region is most likely pop- 
ulated by C stars (iGuUieuszik et al.l l2008l : iMenzies et al.l 
I2OO8I : IWhitelock et al.ll2009l ). although in the Milky Way 
and the Magellanic Clouds it may be populated by other 
objects, such as M-type (O-rich) M i ra variable stars and 
OH/ I R stars (iBessell fc BrettI Il988l : iNikolaev fc Weinberg! 
I2OOOI : iMenzies et al.ll2002l '). Mid-infrared spectroscopy with 
Spitzer Space Telescope of AGB stars in Local Group galax- 
ies has shown that at low metallicity the dusty winds and 
mass-loss rates of C stars are similar to those of Galac- 
tic AGB stars, while mass loss is 1-2 orders of magnitude 
lower in O-rich AGB stars (see iLagadec fc Ziilstral |2008| . 
for a discussion). For example, this is the case for For- 
nax dSph, which has a metallicity similar to that of Leo I 
l|Matsuura et al.ll2007l : iLagadec et al.ll2008l ). All previously 
identified C stars in Leo I are indeed found in region 4. 
NIR photometry for the candidate C stars selected from 
our two-colour diagram is given in Table [5] along with their 
cross-identification with previ ous surveys (unfortunately , no 
matching was possible with lAaronson fc Mould! [l985l be- 
cause finding charts or coordinates are not provided). Most 
of them are previously identified C stars, while 6 are new 
candidates. Note that star No. 11242 (ALW-17), which is 
marginally included in region 4, is not a C star according 



Figure 8. Selection of a clean sample of RGB stars in Leo I using 
the NIR two-colour diagram. In the left panel, we plot stars in 
the upper 1.3 mag of the RGB of Leo I dSph (16.2 < K < 17.5). 
The same magnitude interval is shown in the right panel for the 
external field. 



to !Azzopardi et al] (!l986! ). A few very red stars in region 4 
were not identified in optical bands. For all stars with miss- 
ing optical photometry, visual checks on our EMMI images 
ruled out possible mismatches and confirmed that the opti- 
cal counterparts fall below the detection threshold. A special 
case is star No. 8336 (also M02-D), which is very bright in 
the optical {V = 18.7) compared to other stars in Table O 
and has a blue optical colour B ~ V — 0.04, inconsistent 
with that of C stars. Inspection of both optical and NIR 
images confirmed its star-like shape. This object is likely to 
represent a chance supe rposition or photometric blend (see 
iGullieuszik et al.l!2007b! . for a similar case in SagDIG). 

Figure [71 shows the location of Leo I stars selected using 
the two-colour diagram in the CMD. The O-rich AGB star 
candidates from Tabled! are distributed on a nearly vertical 
sequence above the TRGB, overlapping in magnitude and 
colour with the faintest C stars. Note that most of the AGB 
stars within 1 mag above the TRGB are O-rich stars. The 
stars in Tabl e [5! closely follow th e mean locus of C stars in 
LG galaxies jTotten et al.!!2000l ) up to J-i<' ~ 1.7. Those 
with redder colours show a large scatter consistent with 
the known variability of upper- AGB c arbon stars and the 
onset of circum stcUar dust envelopes jMenzies et al.![2002! : 
!Whitelock fc Fe ast 2000J. 

We note incidentally that the same two-colour tech- 
nique can be used for discriminating RGB stars in Leo I 
from the Milky Way foreground. This is shown in Fig. |8l 
where we plot the NIR two-colour diagram of stars in the 
magnitude interval 16.2 < K < 17.5, which corresponds 
to the upper 1.3 mag of the RGB of Leo I. The diagram is 
plotted for both the Leo I field and the control field, and is 
slightly more scattered than the data in Fig.[6lbecause of the 
increased photometric error. The stars populating region 3 
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Figure 9. Upper panel: comparison of tlie cumulative radial 
distributions of RGB stars and the samples of C- and O-rich AGB 
stars, as defined by our selection criteria based on the NIR two- 
colour diagram. Lower panel: radial variation in the number of 
AGB stars, normalised to the number of stars in the upper 2 mag 
of the RGB. Stars were counted in elliptical annuli with a Of4 
step in semi-major axis (see text for details). The error bars were 
obtained by error propagation of Poisson errors on stars counts. 



Figure 10. Change in the AGB star population as a function 
of the distance from the centre of Leo I. Left: CMD of stars in an 
inner elliptical region of Leo I with semi-major axis r = 2.' 5. A box 
shows our selection of RGB stars. O-rich AGB stars (indicated by 
the letter M, squares) and C stars {circles) were selected from the 
NIR two-colour diagram. Right: The same, except for an elliptical 
ring between r = 2.'5 and r = 6.'0. Outside the central region, the 
fraction of C stars relative to the RGB stars is lower by a factor 
of 3. 



are mostly Galactic dwarfs. This is confirmed hy a. {V — K), 
K plot of stars in the field of Leo I, where the wide baseline 
of the combined optical-NIR colour allows a neat separation 
of the foreground stars. 

The great majority of stars in regions 1 and 2 are red 
giants in Leo I. The near absence of stars in region 1 for 
the comparison field is easily explained by the fact that 
giant stars in the Milky Way are much brighter than the 
magnitude range of the RGB of LeoL While the stars in 
region 1 are mostly in Leo I, the region 2 sample can be 
shghtly contaminated by distant Galactic dwarf stars (but 
in the case of Leo I the contamination is negligible) . The two 
samples of RGB stars should probably be grouped together 
to avoid a bias in colour and metallicity. In conclusion, a 
selection based on broad-band JHK imaging alone proves 
to be a valuable tool for selecting red giants stars. The in- 
creasing availability of wide-field NIR imagers makes this 
technique an interesti ng complement to int ermediate-band 
optical searches (e.g.. [Maiewski et al.ll2000l ) to define clean 
samples of member stars in resolved systems for follow-up 
spectroscopic studies. 

6.2 Radial gradient in the AGB population 

Leo I is one of the few dSph for which no evidence of a pop- 
ulation gradient was known from the literature. Recently , 
our spectroscopic study of Leo I jGuIlieuszik et al] |2009| ) 
gave a marginal detection of a radial metallicity gradient 
among Leo I RGB stars, in contrast with the conclusions of 
iKoch et al.l (|2007l ). We also found evidence for a concentra- 



tion of younger stars in the central region of Leo I. Then it is 
important to confirm the presence of stellar population gra- 
dients using AGB stars as tracers of intermediate-age stellar 
populations. Radial distances were therefore calculated for 
AGB and RGB stars assuming that the centre of Leo I is lo- 
cated at a = 10''08™26=68, 5 = -f 12°18'19'.'7 (J2000) and the 
stellar density distribution is approximately ell iptical with 
ellipt icity e = 0.37 and position angle P.A.=84° (|Sohn et al.l 
[2003) • The radial distribution of AGB stars, normalised to 
the number of stars in the upper 2 mag of the RGB, is plot- 
ted in Fig. |9] (lower panel). C- and O-rich AGB stars were 
counted in elliptical annuli with a 0'4 step in semi-major 
axis. The errors on the relative number of AGB stars were 
calculated assuming Poisson statistic for the star counts. A 
central concentration of AGB stars relative to RGB stars is 
apparent in the plot, with a decline of a factor of ~ 3 at 6' 
with respect to the centre. The statistical significance of this 
result is larger when we consider C stars only (Fig. |9] upper 
panel). A Kolmogorov-Smirnov (KS) test used to compare 
the cumulative distributions of AGB stars (C- and O-rich) 
to that of RGB stars gave the following results. The null hy- 
pothesis that AGB and RGB stars are drawn from the same 
parent population can be rejected with a 95% confidence; 
the probability rises to 97% if we consider only C stars. For 
the O-rich AGB stars, the null hypothesis cannot be rejected 
at any significant level. These tests imply that the C-type 
AGB stars show a statistically significant central concentra- 
tion, whereas the spatial distribution of O-rich AGB stars is 
similar to that of the bulk population of red giant stars. 
The central excess of C stars is also evident in Fig. 1101 
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where we compare AGB stars in an inner 2f5 region of Leo I 
and in an elliptical ring with 2^5 < r < 6'.0. While the frac- 
tion of O-rich AGB stars relative to the number of RGB 
stars is constant in the two regions, the fraction of C-stars 
is 3 times higher in the central region. Again, we conclude 
that the difference in the radial distribution of AGB and 
RGB stars is driven by an excess of C stars, which trace 
the intermediate-age stellar population in the age range 1-5 
Gyr. This confirms the presence of an age gradient in Leo I 
as su ggested by spectroscopy of RGB stars l|Gullieuszik et ahl 
I2OO9I ). Indeed, the number of C-stars is much more depen- 
dent on age and SFH than the number of O-rich AGB and 
RGB stars, as illustrated in Fig. [11] and discussed in more 
detail in the next section. Also, the mean age of red clump 
stars in Leo I is ~ 5 Gyr, definitely older than that of lu- 
minous AG B stars, whi c h exp lains why no age gradient was 
detected bv lHeld et all (|2000l ). 



7 COMPARISON WITH THEORETICAL 
MODELS 

The present observations sample quite completely the AGB 
population of Leo I in the surveyed area, except perhaps for 
th e stars most absorbed by circumstellar dust. As discussed 
bv lGuUieuszik et al.l (|2008l ) for Leo II, this kind of data is im- 
portant for testing thermally-pulsing AGB (TP- AGB) mod- 
els in the regime of low metallicity. In the case of Leo I, the 
presence of recent star formation makes it a particularly in- 
teresting case for a study of AGB stars of different ages and 
masses. 

Since the HST observations used in the determination 
of the SFH of Leo I were obtained in a central region, and 
given the presence of a significant gradient in the stellar 
populations of Leo I, we performed our analysis on the in- 
ner (projected) 2.5 arcmin of the galaxy. Inside this area, 
the A^AGB /N-R.GB ratio is almost constant and comprised be- 
tween 0.10 and 0.12. Selecting this small region also dramat- 
ically reduces the number of foreground stars, whereas the 
number of Leo I stars is reduced by just a factor of 2 relative 
to the total observed area (see Fig. llOf) . 



7.1 Simulating the photometry 

The Leo I star counts can be comp ared with the recent set o f 
TP- AGB evolutionary tracks from lMarigo fc Girardil (|20o3). 
The procedure followed t o simulate the Leo I obse rvations is 
essentially the same as in lGuUieuszik et all (|2008l '). to which 
the reader is referred for details. Suffice it to mention that 
the l atest version of the TRlLEGAl[f| population synthesis 
code iGirardi et al.ll2005l : fGirardi fc Marigdl2007h was used 
to simulate both the Milky Way foreground and the Leo I 
dSph. For this simulation, we assumed for Leo I a distance of 
250 Kpc and zero reddening in the near-infrared. The sim- 
ulations of the TP- AGB include important effects such as 
luminosity and temperature variations during thermal pulse 
cycles, and the p rocessing of radiation by circumstellar dust 
according t o the Groenewegen (|2006l ) prescriptions (details 
are given in iMarigo et al.ll20oi ) . 
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Top panel: The SF Hs used in this work: 

llDolphin et al.ll2005l . dotted lines) and dCallart et al.lll999l . solid 
lines). Bottom panel: The age distribution (number of stars per 
age interval) of different types of stars as a function of age, for a 
galaxy model forming stars at a constant rate from to 15 Gyr, 
at a constant Z = 0.0005 metallicity. The stellar kinds plotted are 
RGB stars within 2 mag of the TRGB (multiplied by 0.4; dotted 
line), and both O-rich (dashed line) and C-rich (solid line) giants 
within 2 mag above the oldest TRGB. 



The total mass of simulated Leo I stars inside the area 
of our observations was chosen so as to reproduce the star 
counts in the upper 2 magnitudes of the RGB, for which our 
observations are complete. The relative star-formation rate 
of Leo I was taken from tw o different sources, as depicted in 
the up per panel of Fig . 1111 
9) and IDolphin et al.l 



|http: //stev. oapd. inaf . It/trll6gal| 



. iGallart et al.l (|l999l . their figure 
(|2005h . In both cases the SFH was 
derived by inversion of a deep CMD from HST observations. 
Both determinations indicate that the star-formation rate of 
Leo I increased in the last few Gyr. This feature appears as 
a broad peak in the SF H between ~ 1 and 7 Gyr in the case 
of iGallart et~aLl |l99i), and a s a much more confin ed peak, 
close to 2.5 Gyr, in the case of lDolphin et al.1 (|2005l ). For the 
Leo I metallicity, we adopted a Gaussian distribution with 
mean [M/H]= —1.3 and dispersion 0.2 dex. This metallicity 
distribution was assumed to be the same for all ages. This is 
a fair approximation since the bulk of the AGB stars arise 
from an age interv al where spectroscopy suggests modest 
chemical evolution (jGullieuszik et aDl2009h . 

To reduce the statistic fiuctuations in the numbers of 
predicted stars, each simulation was run 100 times with dif- 
ferent random seeds. One example of such simulations is 
shown in Fig. 1121 which reproduces quite well the observed 
diagram shown in Fig. [T] The main stellar features are ac- 
counted for by the model, including the colour distribution 
of C stars which, with the new prescriptions for circum- 
stellar dust extinction, turns out to be consistent with the 
extended AGB red tail in the data. The simulations do not 
contain the objects found in the bottom-right part of the 
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Figure 12. One of t he many simulated CMDs for Leo I, using the 
iDolphin et"all bOOSi) SFH. In the electronic version of this paper, 
different colours mark different kinds of stars, namely: Milky Way 
disk (green) and halo (magenta), Leo I pre-TP-AGB stars (dark), 
and Leo I TP- AGB stars both O-rich (blue) and C-rich (red). 

diagram {K > 18, J — K > 1.2 ), which likely correspond t o 
background galaxies (see, e.g., iNikolaev fc Weinberg|[200ol ) . 

7.2 Comparing foreground and RGB counts 

According to the model counts, the expected number of fore- 
ground Milky Way stars in our 3.4 X 10"^ deg^ area (i.e. the 
area of an ellipse with semi-major axis r = 2'.5 and ellipticity 
0.37), in the 13 < A" < 16 magnitude interval, is 2.9 ± 1.7; 
this agrees with the 3 objects observed in regions 2+3 at 
K < 16 and r < 2'.5. In view of this agreement, our simu- 
lations can be used to infer the field contamination in other 
CMD regions. We find that a total of 6.5 ± 2.3 foreground 
stars are expected to contaminate the uppermost 2 mag of 
the RGB in Leo I, in the 16.2 < K < 18.2 interval. The ob- 
served number of stars in this CMD region is 276 (Fig. llOp . 
which largely outnumbers the expected foreground stars. We 
thus conclude that about 270 genuine RGB stars in Leo I are 
expected within 2 mag of the tip. The Leo I simulations are 
then scaled so as to obtain a total of 270 ± 18 RGB stars in 
the upper 2 magnitudes of the RGB. 

7.3 Comparing AGB counts 

Using 100 simulations for each SFH, we obtained the 
following resu l ts for AGB star counts in Leo I. For the 
iDolphin et al] (|2005h SFH, we expect to find 14.8 ± 4.1 O- 
rich giants above the RGB tip. Among them, 11.6 ± 3.7 are 
genuine TP-AGB stars, the remaining being either early- 
AGB or core helium burning stars of the youngest ages. The 
predicted C-rich AGB stars are 40.6 ±6.9 in total, 30.2 ±5.9 
brighter than the TRGB in the K band. In the inner 2!5, the 
data present 11 O-rich and 17 C-rich stars in regions 1 and 
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Figure 13. Simulated LPs for luminous stars in Leo I {solid 
line), separated as O-rich giants above the TRG B (top panel) 
and C -rich TP-AGB stars (bottom panel), for the iGallart et aU 
lll999l) SFH. The shaded histograms correspond to all candidate 
AGB stars in Tables |4] and \5\ Also shown for comparison is 
the simulated LP for the SMC, both including (dashed line) and 
excluding {dotted line) stars younger than 1.5 Gyr. 



4 of Fig. |6l While the number of O-giants in the simulations 
agrees with the observed one, the excess of C-rich giants is 
significan t, by a fac tor of about 2. Alternatively, using the 
ICallart ct all (|l999t ) SFH, we predict 10.5 ± 3.4 O-rich gi- 
ants above the RGB tip (7.0±2.8 genuine O-rich TP-AGB), 
and 47.2 ± 6.9 carbon stars (34.0 ± 5.7 above the TRGB). 
These numbers indicate a small deficit of O-rich stars (not 
significant, however), and again an excess in the predicted 
number of C stars by a factor of ~ 2. 

The mean C-star luminosity function is compared to 
the data in Fig. [13] (bottom panel) . As suggested by a com- 
parison between the data in Fig. [T] and the simulations in 
Fig. 1121 the simulated C stars are fainter on the mean than 
the observed stars. A KS test indicates a negligible probabil- 
ity that the observed distribution is drawn from the model 
distribution. In particular, we note that the predicted num- 
ber of C stars fainter than the RGB tip, which is significant 
in our models, contrasts with the paucity of faint C stars in 
the spectroscopic surveys (although the completeness of the 
observational searches of C stars for objects fainter than the 
RGB tip is not well understood). 

A similar comparison for O-rich AGB stars above the 
RGB tip (Fig. 1131 top panel) shows that the models pre- 
dict a distribution extended to brighter magnitudes than ob- 
served, with a KS probability that the observed distribution 
is drawn from the theoretical model of just 1.5%. Although 
the number of O-rich AGB stars seems to be correctly pre- 
dicted, the stars are brighter than observed by ~ 0.5 mag. 

Thus, as in the case of Leo II, there are significant 
discrepancies between the LFs predicted by current low- 
metallicity TP-AGB models and the observed data. How- 
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ever, the trends are quite different in the two cases. In Leo II 
the lifetimes of O-rich AGB stars are overestimated, while 
the C stars are rep roduced in seemingly acceptable numbers 
IIGullieuszik et al. I [2008). For Leo I, the situation regarding 
O-rich and C-rich AGB star counts is the opposite. The dif- 
ference is likely related to the different mean ages of AGB 
stars in these two galaxies, since the Leo II stellar popula- 
tions are significantly older on average than those in Leo I. 
The few observed C stars in Leo II sample the lower limit of 
the mass interval for which C stars can be formed by thermal 
pulses, with masses close to 1.0 Mq and ages between ~ 5 
and 7 Gyr. Conversely, the numerous C stars in Leo I sam- 
ple the M > 1.3 Mq interval for ages younger than 3 Gyr 
[Z = 0.001). The differences found for C stars are therefore 
not surprising. 

7.4 Comparison with the SMC 

The current models of TP- AGB stars were calibrated using 
constraints from the Large and Small Magellanic Clouds, in 
particular their global C-star luminosity f unctions, and AGB 
lifetim es derived from their star clusters (jMarigo fc GirardU 
Therefore, the behaviour of these models reflects 
constraints set by stellar populations with metallicity be- 
tween [Fe/H] = —1.0 and —0.4. Leo I has a mean metallicity 
~ —1.3 dex, which is 0.5 dex lower than that of the young 
SMC populations, and comparable to that of intermediate- 
age and old populations in the SMC. Thus, it is not clear 
why models calibrated on the SMC observations predict the 
numbers and luminosities of the Leo I AGB stars so poorly. 
Is the moderate difference in metallicity the main reason of 
the observed discrepancies ? 

To clarify these questions, we made two tests. First, 
we tested the effects of changing the metallicity on the LF. 
Simulations were produced with metallicity varying from 
[M/H] = -1.6 to -1.0. The resulting LFs were nearly the 
same as those obtained for [M/H] = —1.3, indicating that 
the problem cannot reside in the metallicities. Second, we 
compared the properties of AGB stars in Leo I and the SMC. 
New SMC simula t ions w ere performed using the SFH of 
[Harris fc Zaritskv (|2004|) together wit h the age- metallicity 
relation of lPagel fc Tautvaisiend (jl998h . We used the trile- 
GAL code in its current version, yielding more detailed simu- 
lation s than the original code presented in lMarigo fc Girardil 
The new simulations, which have the same features 
(e.g., detailed modelling of the photometry) and the same 
physics as used for Leo I, essentially confirm the results of 
iMarigo fc Girardil l|200'if l concerning the shape of the C-star 
LF and the good overall fit to the SMC data. This SMC 
model is shown in Fig. [13] (dashed histogram), shifted to 
the distance and reddening of Leo I dSph. The C stars in 
the SMC are apparently much brighter on the mean than 
those in Leo I (continuous histogram, using the Gallart et 
al. SFH). The dotted histogram shows the same model for 
the SMC, but excluding the stars with ages younger than 

1.5 Gyr. It is then clear that the bulk of the bright C stars 
in the SMC comes from the stron g star format ion episode 
experienced in the recent past (see [Harris fc Zaritskv 20o3) . 

The C-star LF is clearly a sensitive function of the SFH 
of the system. This is confirmed by Fig. 1141 where we present 
a colour-coded plot of the number of AGB stars in Mk and 
age bins in the case of a constant star-formation rate and 
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Figure 14. Simulated number density of evolved (RGB and 
AGB) star in a stellar system with constant star-formation rate 
and metallicity [M/H] = —1.3. The number of stars is colour- 
coded as a function of age and ii'-band absolute magnitude. Top 
panel: stars that did not undergo any thermal pulses. These 
stars comprise RGB stars and early- AGB stars (older than ~ 1 
Gyr) and red supergiants evolving from intermediate-mass stars 
younger than 1 Gyr. Middle panel: luminosity distribution of 
C-rich AGB stars. Bottom panel: O-rich TP-AGB stars. 



fixed metallicity equal to that of Leo I. In Fig. 1141 most of 
the TP-AGB stars younger than ~ 6.5 Gyr become C stars 
soon after they enter the TP-AGB phase, except for some 
TP-AGB stars in the 2-5 Mq range which also appear as O- 
rich because either are evolving through the first part of their 
TP-AGB phase before being converted to C stars, or (those 
younger than 0.5 Gyr) are massive enough to experience 
hot-bottom burning, thu s keeping their surface C/0 ratio 
< 1 (e.g., iMarigol 120071 ). As seen in the middle panel of 
Fig. 1141 a sustained star formation activity in the last 1-2 
Gyr does imply a high production rate of C stars. Since a 
recent episode of star formation (^1 Gyr) is absent in Leo I 
dSph, very luminous C stars are missing in the simulations. 

While the C stars in the SMC mostly sample the proper- 
ties of young C stars of metallicity close to —0.4, with initial 
masses higher than ~ 2 Mq , the AGB stars in Leo I mainly 
arise from intermediate-age stars with metallicity —1.3 and 
initial masses close to 1.5 Mq. 



7.5 Models vs. observations 

Then, for the bright end of the C-star LF, it is conceivable 
that the star-formation histories used to calculate the model 
LF are biased towards stellar populations which are on av- 
erage too old. However, it appears likely that the models 
actually overestimate the lifetime of low- and intermediate- 
mass C-type stars, and/or underestimate their luminosities. 
There are different ways by which the lifetimes/luminosities 
may be changed: by increasing the mass-loss rates for C 
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stars, by decreasing the efficiency of dredge-up events (so 
as to more rapidly increase the core mass and luminosity), 
or by efficient dimming of their near-infrared light by cir- 
cumsteUar dust sheUs. Another aspect that might be po- 
tentially relevant to this issue is the fact that at the low 
metallicities considered here the TP-AGB models attain a 
surface C/O > 1 already during the first few sub-luminous 
thermal pulses, before reaching the core-mass luminosity re- 
lation characteristic of the fully-dev eloped TP-AGB regime 
IWagenhuber fc GroenewegenI 19981 ) . Though the effect of a 
very efficient third dredge-up on the stellar luminosity dur- 
ing these earlier stages is still not well assessed, we find 
in the literature some theoretical indication that TP-AGB 
stars might be brighter than commonly assumed by syn - 
thetic models (see, e.g., iHerwig et cj]|l998l : lMowlavilll999t ). 
While exploring the possible solutions requires the calcula- 
tion of new grids of TP-AGB evolutionary models, which is 
beyond the s cope of this paper, the present results (along 
with those of iGuUieuszik et al.ll20oj ) demonstrate that, de- 
spite the possible uncertainties in SFHs and age-metallicity 
relations in nearby dwarf galaxies, the observed properties 
of AGB stars in dSph galaxies provide valid constraints to 
AGB models at low metallicities. 



8 SUMMARY AND CONCLUSIONS 

We have presented new infrared JHK photometry of Leo I 
dSph obtained with the wide-field imager WFCAM at the 
UKIRT. Our NIR photometry, combined with optical B and 
V band photometry, provides information on the spectral 
energy distribution of AGB and RGB stars in Leo I on a 
wide colour baseline. 

We provide a catalogue of JHK PSF-fitting photom- 
etry of evolved stars in Leo I down to K ^ 20.5, over 
a 13!5 X 13!5 area, from which the basic parameters of 
Leo I, distance and metallicity, were re-derived. We mea- 
sured the distance to Leo I from the J, H and K magni- 
tude of the TRG B, using a calibrati on based on Galactic 
globular clusters jValenti et al.|[2[)04br i. Taking into account 
the SFH of Leo I by mean of synthetic CMDs, we corrected 
our measurements for the age difference between Leo I and 
the template GCs. This yielded a corrected distance mod- 
ulus (m — M)o = 22.04 ±0.11, in agreement with previous 
optical results. The V — K colours of RGB stars were used 
to infer a metallicity distribution. We obtained a mean age- 
corrected metallicity [M/H] = -1.24±0.05 (int)±0.15 (syst), 
consistent with recen t results from spe c troscopic analysis 
of Leo I RGB stars (iBosler et al.l l2007l : iKoch et al.l l2007l : 
iGullieuszik et al I bood ). Astar-bv-star comparison of metal- 
licities derived from spectroscopy and V — K colours gave 
residuals generally consistent with th e ages of the s tars esti- 
mated from optical colours (Gullicu szik et al.ll2009l ). There- 
fore, we argued that NIR photometry allows an accurate de- 
termination of the main parameters (distance, metallicity) 
of resolved galaxies provided that a guess of their SFH can 
be made. This conclusion is important for deriving distances 
and metallicities for resolved galaxies out to the Virgo clus- 
ter of galaxies with the future generation of NIR instruments 
(JWST and adaptive optics at Extremely Large Telescopes). 

Using our NIR photometry, we obtained a nearly com- 
plete (with the possible exception of heavily obscured stars) , 



homogeneous sample of AGB star in Leo I, confirming the 
presence of some very red objects not detected in opti- 
cal images, probably luminous thermally-pulsing AGB stars 
reddened by dust e nvelopes (4 out of 5 in common with 
iMenzies etal] |2002| '|. We showed that, if the photometric 
data are of high precision and accuracy, the NIR two-colour 
diagram is able to discriminate between 0-rich and C-rich 
AGB stars brighter than the RGB tip, and against the fore- 
ground Galactic dwarf stars. Our selection criteria for C 
stars brighter than the TRGB are entirely based on NIR 
photometry, hence not affected by the various biases typi- 
cal of optical methods. The same method can be employed 
to select a clean sample of red giant stars by significantly 
reducing the contamination by Galactic dwarf stars. As a 
result, we found 37 C-stars in Leo I, 6 of which are new can- 
didates. We were able to define, for the first time, a sample 
of 39 bona-fide 0-rich AGB stars against the contaminating 
Milky Way foreground, on the basis of their NIR colours. 
Our selection was tested (and confirmed) on a control field. 
From these samples, we derived the luminosity functions of 
O-type and C-type AGB stars in Leo I. 

The AGB stars in our Leo I sample were used as tracers 
of their parent intermediate-age stellar populations. We re- 
vealed a radial gradient in the number of C-type AGB stars 
relative to RGB stars, which implies an increasing fraction 
of intermediate-age (1-3 Gyr) stars in the inner regions of 
the galaxy. The fraction of C stars within an ellipse with 
semi-major axis r = 2f5 (186 pc), having ellipticity and po- 
sition angle appropriate for Leo I, is larger by a factor of 
3 than in the outer field. This confirms the presence of an 
age gradient in the stellar p opulations of Leo I, as s uggested 
by our spectroscopic study (|Gullieuszik et al]|2009l ) . In con- 
trast, the O-rich AGB stars follow the profile of RGB stars. 

The NIR photometry of AGB stars in Leo I was com- 
pared with the prediction of theoretical models. As in the 
case of Leo II, we found significant discrepancies between 
the present TP-AGB models with low metallicity and the 
observed data, with the simulated LFs for C- and 0-rich 
AGB stars not reproducing the observations. In particular, 
the predicted number of C stars fainter than the RGB tip 
is larger in our models than found by spectroscopic surveys. 
These discrepancies can be explained by uncertainties in the 
star formation histories used to calculate the model LFs or, 
more likely, by an overestimate of the lifetime of low- and 
intermediate-mass C-type stars, and/or an underestimate of 
their luminosities, in the current models. 
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Table 4. Candidate O-rich AGB stars in Leo I, as selected from the NIR two-colour diagram 
(Fig. |6l region 1). 
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